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Abstract. Ultrarelativistic electron-positron plasmas can be produced in high-intensity laser fields and 
play a role in various astrophysical situations. Their properties can be calculated using QED at finite 
temperature. Here we will use perturbative QED at finite temperature for calculating various important 
properties, such as the equation of state, dispersion relations of collective plasma modes of photons and 
electrons, Debye screening, damping rates, mean free paths, collision times, transport coefficients, and 
particle production rates, of ultrarelativistic electron-positron plasmas. In particular, we will focus on 
electron-positron plasmas produced with ultra-strong lasers. 

PACS. 52.27.Ep Electron-positron plasmas - ll.10.Wx Finite-temperature field theory 
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1 Introduction 

Plasmas, i.e. (partly) ionized gases, are considered to be 
the fourth state of matter after the solid, liquid, and gaseous 
states. The plasma state dominates in the observable Uni- 
verse: 99% of the visible matter is in the plasma state, 
namely in the form of stars and hot interstellar and inter- 
galactic gases. Plasmas emit light due to the excitation of 
atoms and ions and recombination. Plasmas can be pro- 
duced by high temperatures, such as in the sun or fusion 
reactors, by electric fields (discharges), as used for illumi- 
nation in neon tubes or in lightening, or by radiation, such 
as in the Crab nebula where the pulsar in the center emits 
synchrotron radiation. 

Plasmas can be classified according to various aspects: 
Relativistic plasmas, e.g. the electron-positron plasma in 
a supernova explosion, are plasmas in which the ther- 
mal energy fc#T of the plasma particles is of the order 
of their rest mass energy mc 2 or larger. Quantum plas- 
mas, e.g. the degenerate electron component in a white 
dwarf, are plasmas in which the thermal de Broglie wave 
length Xb — h/{mv t h) is of the order of the interparticle 
distance d or larger. Here Vth is the thermal velocity of the 
particles. Strongly coupled plasmas, e.g. the ion compo- 
nent in white dwarfs, are plasmas in which the interaction 
energy between the particles is larger than their thermal 
energy. In non-relativistic plasmas this corresponds to a 
Coulomb coupling parameter Zc = Q 2 /(dkT) > 1. 

Plasmas in Nature, like in comets, in aurorae, in the 
corona of the sun, in lightening, in flames, or in the sun, 
and in technology, like in discharges or fusion reactors, 
cover a wide range of pressures and temperatures. All of 
these plasmas are non-relativistic, classical, and weakly 
coupled systems. 

Electron-positron plasmas (EPPs) are created in the 
presence of strong electric or magnetic fields or extremely 



high temperatures, where massive pair production sets in. 
For example in a supernova explosion temperatures up to 
3 x 10 11 K corresponding to k B T ~ 30 MeV > 2m e c 2 = 
1.022 MeV will lead to an ultrarelativistic EPP. Also in 
the vicinity of magnctars, i.e. neutron stars with magnetic 
fields B > 10 14 G, and in accretion disks around black 
holes EPPs show up. 

Recently the possibility to create ultrarelativistic EPPs 
with high-intensity lasers (I > 10 18 W/cm 2 ) have been 
discussed. For example two opposite laser pulses hitting a 
thin gold foil will heat up the electrons in the foil up to 
several MeV leading to pair creation pQ . 

In the following we will discuss the properties of an 
ultrarelativistic EPP using quantum field theory (QED) 
at finite temperature. We will follow closely the review 
article [5] . We will not discuss the production mechanism 
and equilibration of the EPP here. 



2 Field theoretic description of an 
electron-positron plasma 

Throughout the paper we will use natural units, i.e. h = 
c = ks = 1, as usual in quantum field theory, in which 
all units are given in powers of MeV. The conversion to 
conventional units can be achieved by 1 = he = 1.97 X 
10~ 13 MeV m from which 1 MeV = 1.60 x 10" 13 J=5.08 x 
10 12 m _1 =1.52 x 10 21 s _1 follows. In these units the elec- 



tron charge e = 0.3 corresponding to a fine structure con- 
stant a = e 2 /(4?r) = 1/137. 



2.1 Equation of state 

We will start with the equation of state of an EPP and 
compute it under the following assumptions: 
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1. ultrarelativistic EPP, i.e. T ^> to, 

2. thermal and chemical equilibrium, 

3. equal electron and positron density, i.e. vanishing 
chemical potential, 

4. ideal gas, i.e. no interactions in the plasma, 

5. infinitely extended, homogeneous, and isotropic EPP. 

We will relax some of these assumptions in the follow- 
ing sections. According to these assumptions the distribu- 
tion function of the electrons and positrons is given by the 
Fermi-Dirac distribution 



n F (p) 



1 



eP/ T + 1 



and of the photons by the Bose- Einstein distribution 

1 



n F (p) 



oV/T _ l ■ 



a) 



(2) 



where the momentum p is identical to the energy E of the 
particles in the ultrarelativistic case. It should be noted 
that the photons are in equilibrium with electrons and 
positrons under the above assumptions, i.e. the system is 
actually an electron-positron-photon gas. 

The particle and energy density can be calculated by 
integrating over the distribution functions. The particle 
number density of the electrons and positrons follows from 
integrating over the Fermi-Dirac distribution as 



9f 



d 3 p 
(2tt) s 



n F (p) 



C(3)T 3 = 0.37 T 3 , (3) 



where g F = 4 is the number of degrees of freedom corre- 
sponding to the electrons and positrons in the two spin 
states. Assuming a temperature of T = 10 MeV, we find 
pf = 370 MeV 3 = 4.9 x 10 40 m" 3 . 

The photon density follows accordingly by integrating 
over the Bose-Einstein distribution with gg — 2 degrees 
of freedom corresponding to the two polarization states as 
P^-(2A 2 )C(3)T 3 = 0.24T 3 . 

The energy density of the electron-positron-photon gas 
is obtained from 



e eq = g F 



d 3 p f d 3 p 
^pn F {p)+g B J j2nf pnB ^ P ' 



lln 2 
~60~ 



(27T) 

T 4 = 1.81 T 4 . 



(4) 



where the photons contribute 36 % to the energy density. 
Here the Boltzmann law, e eq ~ T 4 , holds also for the 
fermions because we neglected their masses. 

For T = 10 MeV we find e e « = 3.8 x 10 29 J m~ a . In a 
volume of 10 12 m 3 (corresponding to the size of a neutron 
star) the total thermal energy of the EPP is 3.8 x 10 41 
J, which corresponds to about 10% of the entire energy 
(without neutrinos) released in a supernova type II ex- 
plosion. In a volume of 1 ^m 3 there is still an energy of 
3.8 x 10 11 J contained. 

The Coulomb coupling parameter of the EPP, which 
is a measure for the non-ideal behavior of a plasma [2j , is 



given by Fc = e /(dT), where d ~ p\ 



-1/3 



2.7 x 10 



-14 



m is the interparticle distance. For T = 10 MeV we find 
r c = 5.3 x 10~ 3 which shows that the EPP is a weakly 
coupled plasma. Therefore the ideal gas results for the 
equation of state derived above are a good approxima- 
tion. After all, interactions in the EPP play an important 
role, for example, for the collective behavior of the plasma 
as discussed in the next section and for equilibration of 
the plasma. Obviously, the interactions can be treated by 
perturbation theory. 



2.2 Collective Phenomena 

Collective effects in a plasma are associated with long- 
range interactions within the plasma. Important examples 
are Debye screening and plasma waves. In non-relativistic 
ion-electron plasmas t 4 plasmas these phenomena can be 
described by classical transport theory (Vlasov equation). 
For example, the electron plasma frequency reads 



' 47re 2 p e 



(5) 



and the Debye screening length due to the electrons in the 
plasma 



A, 



47re 2 /? e 



(6) 



where p e is the electron number density, T e the tempera- 
ture of the electron component, and m e the electron mass. 
In an ultrarelativistic plasma with T> ro the masses can 
be neglected and the important scales are the tempera- 
ture T, called the hard scale, and the soft scale eT, which 
determines the collective phenomena as we will see below. 

Interactions between relativistic electrons and positrons 
can be treated by using perturbative QED. This corre- 
sponds to an expansion in the fine structure constant a. 
Most conveniently Feynman diagrams are considered from 
which via Feynman rules quantities such as scattering 
cross sections, decay and production rates, or life times 
can be calculated directly. In an EPP the interactions take 
place in the presence of a heat bath. Hence we have to 
consider QED at finite temperature. For this purpose the 
Feynman rules are extended to finite temperature which 
can be achieved by using the imaginary or real time for- 
malism [5j[6] . The calculations are similar to the ones done 
already within the last 30 years for the properties of the 
quark-gluon plasma using perturbative QCD at finite tem- 
peratures. As a matter of fact, many results from the 
quark-gluon plasma (see e.g. Ref.JT]) can be directly car- 
ried over to the EPP. 

An important quantity is the polarization tensor or 
photon self-energy. The lowest order diagram for the po- 
larization tensor is shown in Fig.l. Assuming the external 
four momentum K — (fco, k) to be soft, i.e. the frequency 
fco = lj and k = |k| to be much smaller than T, and the 
internal loop momenta to be hard, an analytic result can 
be found using the real or imaginary time formalism 



n L (w,k) = -3r 



'ph 



cj uj + k 
1 -In - 

2k uj — k 
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Fig. 1. One- loop polarization tensor. 
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(7) 



where m p h = eT/3 is called the effective photon mass. For 
T = 10 MeV we get m ph = 1 MeV. 

The crucial quantity from which the collective phe- 
nomena are derived is the dielectric tensor relating the 
macroscopic electric field Di in the medium to the exter- 
nal field Ei (i = x, y, z), i.e. in momentum space 



A(w,k) = 5^Q i (w,k)£; i (w,k). 



(8) 



In the case of an isotropic medium it depends only on lo 
and on k and has two independent components 



dij(w,k) = e T {u,k) 8, 



k% kj 



+ e L (uj,k) 



k^ kj 



(9) 



The dielectric tensor is closely related to the polariza- 
tion tensor or photon self-energy by [10] 



e L (u, k) = 1 - 
e T (u>, k) = 1 - 



n L {uj,k) 
k 2 ' 

n T {m,k) 



where II l and Ut are the longitudinal and transverse 
components of the polarization tensor, respectively. 

The dielectric functions following from (| 10[) and J?} 
can also be derived from the classical Vlasov equation to- 
gether with the Maxwell equations [TT], since the high- 
temperature limit corresponds to the classical limit. 

The dispersion relations of collective plasma modes, 
i.e. propagation of electromagnetic waves in the plasma, 
can be found by using the Maxwell equation, leading to 



e L (io, k) = 0, 
k 2 



(11) 



Combining (JTUJ), (0, and (fTTj) gives the dispersion re- 
lations UL.T(k) of the transverse as well as longitudinal 
plasma waves as shown in Fig. 2. The longitudinal branch, 
which does not exist in vacuum, is called plasmon as in 
the case of non-relativistic plasmas. The transverse branch 
does not play a role in non-relativistic plasmas but is 
equally important as the longitudinal one in relativistic 
plasmas. Both branches start at the plasma frequency 




Fig. 2. Photon dispersion relation. 



ujpi — uiL^ik = 0) = m p h- Consequently the collective 
plasma waves have soft momenta of the order eT. At high 
momenta k 3> m p h the transverse mode approaches the 
free dispersion lot — k, corresponding to a real photon in 
vacuum, whereas the longitudinal mode disappears, i.e. its 
spectral strength is exponentially suppressed. For T = 10 
MeV we find oo pl = 1.5 x 10 21 s" 1 . 

Another important quantity which can be derived from 
the polarization or dielectric tensor is the Debye screen- 
ing length, entering the Yukawa potential of a heavy, non- 
relativistic test charge in the EPP. The Debye screening 
length is given by the static limit of the longitudinal com- 
ponent of the polarization tensor l/iTj,(w = 0), leading 
(10) to \ D = l/( v / 3m p/l ), which is 1.1 x 10~ 13 m for T = 10 
MeV. 

Finally from (J7J) we see that the polarization tensor 
and the dielectric function become imaginary for co 2 < k 2 , 
i.e. below the light cone lo — k, corresponding to Landau 
damping. We also observe that the plasma waves calcu- 
lated at lowest order perturbation theory are undamped 
since they are located at lo > k. 

A complete new phenomenon that does not appear 
in non-relativistic plasmas is the existence of fcrmionic 
plasma waves because all fermion masses are much too 
large in the non-relativistic case. Their dispersion relations 
follow from the pole of the electron propagator containing 
the electron self-energy. Using again the high tempera- 
ture approximation for the one-loop electron self-energy 
of Fig.3 leads to (P = (p , p), p = |p|) [8E2] 



2{P) = -o(po,p) P"7m ~ &(Po,P) 7o 



(12) 



with 



a{Po,p) 
b{Po,p) 



1 

br 
1 

Ap 2 



[tr(P^E)-p tr( l0 E)} , 



[P 2 tr( 10 Z)-Potr(P^E)] , (13) 
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Fig. 3. One-loop electron self-energy. 

where the traces over the 7 matrices are given by 
tr(P^ fl £)=4m F , 

tr(j S)=2 ml - In ^_±P (14) 
P Po-P 

with the effective electron mass niF — eT /■*/&, which is 
1.1 MeV at T= 10 MeV. 

The full electron propagator in the hclicity represen- 
tation is given by |13j 



S*(P) 



where 



1 



2D+(P) 



(lo-p-7) + 



1 



2D_(P) 



(lo+p-l) , (15) 



3 




p/m F 

Fig. 4. Electron dispersion relation. 



D±(P) = -Po±P+^ [±tr(P^ E) ~ (±P0~P) M70 S)] 

(16) 

The dispersion relations following from the pole of this 
propagator are shown in Fig. 4. Again two branches show 
up, one with a positive ratio of the helicity to chirality 
(w+) following from D + = 0, the other one with a neg- 
ative ratio (w~) following from _D_ = 0, called plasmino 
[13 . The plasmino branch uj_, which does not exist in vac- 
uum, shows an interesting behavior, namely a minimum at 
k = OAlrriF, which may lead to van Hove singularities [13l 
[T4] . Whether these van Hove singularities will lead to ob- 
servable effects in the EPP, e.g. in the electron spectrum, 
is a very interesting question which should be investigated 
in detail. It could open the exciting possibility to observe a 
new collective plasma wave, the plasmino, experimentally 
in a laser induced EPP. 



2.3 Transport properties 

Now we want to consider the interaction and properties 
of particles in the plasma with hard momenta, i.e. of the 
order of T or larger. In particular we are interested in 
damping and transport rates, mean free paths, collision 
times, energy losses of these particles and other transport 
properties such as the shear viscosity of the EPP. 

It was shown by Braaten and Pisarski [15] that a con- 
sistent treatment of gauge theories such as QED at finite 
temperature, i.e. for obtaining results that are gauge in- 
dependent, infrared finite, and complete to leading order, 
require the use of an effective perturbation theory using 
resummed Green functions based on the hard thermal loop 
(HTL) approximation (HTL resummation technique). 



The damping rate of an electron or positron in the 
EPP is defined as the imaginary part of the dispersion re- 
lation ujl.t{p)- To lowest order it follows from the elastic 
scattering diagram of Fig. 5. In the case of a hard electron 
or positron with momenta of the order of T or higher it 
exhibits a quadratic infrared (IR) divergence which can be 
reduced to a logarithmic one using a HTL resummed pho- 
ton propagator. This logarithmic singularity is expected 
to be cut-off by higher order contributions leading to [7] 



7e 



e 2 T 1 
— — ln- 

4-7T e 



(17) 



within logarithmic accuracy, i.e. the constant under the 
logarithm is not determined. For T = 10 MeV we obtain 
7 e = 86 keV, which is much smaller than lo p i — 1 MeV, 
showing that the EPP is not overdamped. 

Physically more important are the transport rates P 
which are related to the mean free path and collision time 
of electrons and positrons in the EPP. They differ from 
the damping rate in cutting off the long range interactions 
with small scattering angles 9 by a factor (1 — cos 8) under 
the integral defining the rate [I] . This leads to an improve- 
ment of the IR behavior (logarithmic instead of quadratic 
singularity in perturbation theory) and a finite result using 
the HTL method. Logarithmic divergent quantities can be 
treated consistently by splitting them into a soft part and 
a hard part, where the soft part is calculated using the 
HTL resummation technique [16] . For the transport rate 
we find to logarithmic accuracy 



e 4 T 3 1 
37rs e ' 



(18) 
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Fig. 6. Diagrams defining the collisional energy loss. 



K 



Fig. 5. Lowest order diagram for electron-electron scattering. 

where the Mandelstam variable s = (P + K) 2 is the square 
of the sum of the four momenta of the incoming particles 
in the scattering diagram of Fig. 5. For thermal particles we 
replace s by its thermal average (s) — 2(p) e (k) e ~ 19. 3T 2 , 
where (p) e — (k) e — e e e q / p e e q = 3.1 IT. Assuming again 
T = 10 MeV, we get P e = 0.54 keV. 

The mean free path A™-^ p and collision time r e of the 
plasma particles (electrons and positrons) are given by the 
inverse of the transport rate 1/-T e , leading to A™^ p = 0.37 
nm and r e = 1.2 x 10~ 18 s at T = 10 MeV. 

In a non-relativistic plasma the shear viscosity can be 
estimated from elementary kinetic theory as |17j 



I ^2 Pi (Pi) K 



mfp 



(19) 



where the sum is performed over the various components 
of the system. In an relativistic plasma the coefficient 1/3 
should be replaced by 4/15 [18 . Using the mean free path 
following from (fT8|) . the density of (|3|), and the thermal 
momentum (p) e — 3.1 IT, the shear viscosity is given by 
(within logarithmic accuracy) 



55.8r 3 
e 4 ln(l/e)' 



(20) 



At T = 10 MeV the shear viscosity coefficient is r\ e = 
7.9 x 10 10 Pa s. 

Another quantity of interest in a plasma is its stopping 
power or the energy loss of an energetic particle in the 
plasma. There are two contributions, namely the energy 
loss by collisions and the radiative one by bremsstrahlung. 
In a relativistic plasma the latter one becomes important. 
The collisional energy loss is given by the mean energy 
transfer divided by the mean free path leading to [19] 



dE 
dx 



dj tu 



(21) 



where v is the particle velocity, 7 the damping or inter- 
action rate proportional to the plasma density and the 
collision cross section, and u the energy transfer from the 
energetic particle to the plasma particle in the collision. 
Using for the collision cross section the lowest order dia- 
grams in Fig. 6, the collisional energy of a muon with mass 
M in an EPP has been calculated by Braaten and Thoma 
|19j applying the HTL resummation technique 



dE 
dx 



e 4 T 2 
24tt 



1 - v 2 



In 



1 



1 - v 



ln J! + ln 7 + A{v] 

(22) 



/\/\/\/\ 



/\/\/\/\ 



Fig. 7. Diagrams denning the photon damping rate. 



where A(v) is a slowly varying function of the muon ve- 
locity v between 1.3 and 1.5. 

The collisional energy loss of an electron with energy 
E 3> T is approximately given by [19] 



dE 
dx 



e 4 T 2 , Ib.ZE 
In ■ 



48tt 



e 2 T 



(23) 



This leads to an energy loss of 200 MeV/nm for an electron 
or positron with an energy of E — 100 MeV at T = 10 
MeV, showing that such an electron is stopped (thermal- 
ized) within a fraction of a nanometer. So far no calcu- 
lations of the radiative energy loss in an EPP have been 
performed to our knowledge. 

The damping rate of a photon in an EPP follows from 
the diagram in Fig. 7, where a HTL resummed electron 
propagator has to be used in case of soft momenta of the 
exchanged electron (positron). In contrast to the electron 
damping rate, the photon rate is infrared finite using the 
HTL method due to the presence of an electron propa- 
gator in Fig. 7 instead of the photon propagator in Fig. 5. 
Hence there is no need to cut off the long range interac- 
tion introducing a transport cross section. The result for 
a photon with energy E — p reads [5D] 



r. 



ph 



e 4 T 2 , 3.88E 
In . 

e 2 T 



64irE 



(24) 



The mean free path and the collision time of photons in 
an EPP are given by l/r p h- For a thermal photon with 

2.75T at T = 10 



the mean momentum (p) p h = t£h 



/Pph ~ 



MeV the mean free path A^/ p = 0.28 nm and the collision 

time T p h = 9.4 x 10~ 19 s follow. Actually the damping rate 
given in (|24|1 is a lower limit as higher order effects will 
enlarge it. As a matter of fact, the photon production rate 
in a QGP, which is the inverse process of the damping rate 
[2"U] . was shown to be about a factor of 2 larger taking 
bremsstrahlung into account [2T] , 

For the viscosity of the photon component using the 
above mean free path, the photon density (see above), and 
the mean photon energy (p) p h = 2.75T we find 



48. 7T 3 



ln(3.27/e) 



(25) 
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corresponding to 3.5 x 10 10 Pa s at T — 10 MeV. Hence the 
viscosity of the EPP 77 = ije+Vph has similar contributions 
from the electrons and photons. 

A more advanced calculation of the total viscosity of 
the EPP based on the Kubo formula yields within loga- 
rithmic accuracy [22] 



V 



188T 3 
4 ln(l/e)' 



(26) 



This result is about a factor of 1.5 larger than the one 
presented here based on the elementary kinetic theory, 
which is typically valid within a factor of 2 [17] . 

A summary of the QED results of the EPP properties 
discussed above is presented in the table below. 



[23j . The fugacity is given by the ratio of the experimen- 
tal to equilibrium particle density, since the experimental 
density follows from integrating over the non-equilibrium 
distribution, i.e. 



9f 



d 3 p 
(2tt) 3 



Xn F (p) = A/?, 



eg 



A = 10" 



(27) 



Using the real time formalism, QED perturbation the- 
ory and the HTL method can also be extended to non- 
equilibrium situations like the one discussed above |24j . 
For example, the effective photon mass is given now by 



l ph 



4e 2 f°° 

3^1 d PP^( P )- 



(28) 



Quantity 


Formula 


Value at 1 — 1U MeV 


Electron-Positron Density 


= 3/7T C(3) T A 
eq _ 9 /_2 /-(o\ t3 
Pph - Z / n ^\ 6 ) 1 

ef = 7tt 2 /60 T 4 
e ;^ 2 /15T 4 
e e 9 = 11tt 2 /60T 4 


4.9 x 10 m 


riioioii ueiisi^y 


O.Z X 1U III 


Electron-Positron Energy Density 


2.4 x 10 29 J m" 3 


Photon Energy Density 


1.4 x 10 29 J m~ 3 


Total Energy Density 


3.8 x 10 29 J m~ 3 


Thermal electron momentum 


(p) e = ef/pf = 3.11T 


31 MeV 


Thermal photon momentum 


(P)ph = e;i/PH = 2.75T 


28 MeV 


Interparticle distance 




2.7 x 10" 14 m 


Coulomb Coupling Parameter 


A = e 2 /{dT) 


5.3 x 10~ 3 


Effective Photon Mass 


m p h = eT/3 


1 MeV 


Plasma Frequency 


U) p i — m p h 


1.5 x 10 21 s" 1 


Debye Screening Length 


X D = l/(V3m ph ) 


1.1 x 10" 13 m 


Effective Electron Mass 


m F = eT/(2%/2) 

7e = e 2 T/(4n) ln(l/e) 


1.1 MeV 


Electron Damping Rate 


86 keV 


Electron Transport Rate 


r e = e 4 T 3 /(37rs) ln(l/e) 

r ph = e 4 T 2 /{MiTE) ln(3.88£;/e 2 r) 


0.54 keV for s = 19.3 T 2 


Photon Damping Rate 


0.70 keV for E = 2.75 T 


Electron Mean Free Path 


X mf P = y Fe 


0.37 nm 


Photon Mean Free Path 


A ph -t/J-ph 


0.28 nm 


Electron Collision Time 


T e = l/r e 


1.2 x 10" 18 s 


Photon Collision Time 


T ph = l/Tph 


9.4 x 10" 19 s 


Electron Viscosity 


f] e = 55.8T 3 /[e 4 ln(l/e)] 
rjph = 48.7T 3 /[e 4 ln(3.27/e)] 


7.9 x 10 10 Pa s 


Photon Viscosity 


3.5 x 10 10 Pas 


Total Viscosity 


T) = r] e + r/ph 


(1.1- 1.6) x 10 11 Pas 


Electron Energy Loss 


dE/dx = e 4 T 2 /(487r) ln(15.3^/e 2 T) 


200 MeV/nm for E = 100 MeV 



2.4 Non-equilibrium and finite chemical potential 

EPPs produced in strong laser fields are probably not in 
complete equilibrium. For example, it has been predicted 
by Shen and Meyer-ter-Vehn p] that a positron density of 
about 5 x 10 28 m -3 at a temperature of 10 MeV can be 
reached. This density deviates from the equilibrium den- 
sity ((3]) by 12 orders of magnitude. In the following we 
will therefore assume that the EPP produced by lasers 
is in thermal but not in chemical equilibrium. Then we 
can replace the distribution functions for the electrons 
and positrons by Fcrmi-Dirac distributions multiplied by 
a fugacity factor A describing the deviation from chem- 
ical equilibrium, Jf{p) = Xxif{p)- This assumption has 
been used for example for describing the chemical equili- 
bration of the QGP in ultrarelativistic heavy-ion collisions 



For T — 10 MeV we then find for the non-equilibrium 
photon mass rn p l< ^ leq = \f\m v h — 1 eV and the plasma 
frequency Lo p l ° neq — 1.5 x 10 15 Hz. The Debye screening 
length in such an EPP is Ad = 0.1 /im. In order to speak of 
a plasma with collective behavior its dimension L should 
be much larger than Ad, i.e. at least of the order of 1 ^m. 

Furthermore, an anisotropic EPP can also be described 
by quantum field theoretic methods [25] . In this case in- 
stabilities can occur [26] , 

Finally a possible difference between the positron den- 
sity and the electron density can be treated by introducing 
a finite chemical potential /i, i.e. using the distribution 

= e (p±M)/r + i ( 29 ) 



Markus H. Thoma: Ultrarelativistic Electron-Positron Plasma 



7 




e |u 

Fig. 8. Lowest order contribution of the muon production. 

for the electrons (negative sign) and positrons (positive 
sign). Such a difference comes from the fact that the laser 
produced EPP is embedded in a hot electron and cold ion 
background of the target. Therefore there will be an excess 
of electrons over positrons in the hot, relativistic EPP. The 
methods described above, such as the HTL resummation, 
can be generalized easily to this case |27j . For example, 
the energy density is given by 

or the effective photon energy by 

2 e 2 T 2 ( 3fi 2 \ 
m 2 ph = — (l + Jjr 2 j . (31) 



3 Conclusions 

The aim of this presentation is the prediction of proper- 
ties of ultrarelativistic thermalized EPPs produced in laser 
fields or supernovae. The ultrarelativistic EPP is a weakly 
coupled system. Therefore its equation of state can be de- 
scribed in first approximation by an ideal ultrarelativistic 
gas. Interactions within the EPP can be described by per- 
turbative QED at finite temperature. In this way collective 
phenomena (plasma waves, Debye screening) and trans- 
port properties (damping rates, mean fee paths, relaxation 
times, production rates, viscosity, energy loss) can be com- 
puted. A complete new phenomenon, fermionic plasma 
waves (plasmino) , which is absent in non-relativistic plas- 
mas, might be observable by van Hove singularities. The 
deviation from chemical equilibrium, as expected for laser 
produced EPPs, can also be treated by extending QED to 
the non-equilibrium case. 
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2.5 Particle Production 

At high temperatures above 10 MeV also other parti- 
cle species will be produced, e.g. muons with a mass of 
rrifj, = 106 MeV. Their rate follows to lowest order from 
the diagram in Fig. 8 (Born term). We assume that m e <C 
T -c wij, holds. The first inequality implies that the elec- 
tron mass can be put to zero and the latter inequality 
implies that muons are not equilibrated. Then the muon 
production rate to lowest order (e~e + — > 7* — > [iT fi + ) is 
given by (for more details see Ref. [2] ) 

dN a 2 f 2m* \ / _ 4m£\ 1/2 T 

d 4 xd 4 p 24tt 4 y M 2 J [ M 2 J p 

- ) , ln l + exp[-(^ + p)/(2T)] 
exp(E/T)-l l + exp[-(£-p)/(2T)]' v ' 

where M 2 = E 2 — p 2 is the invariant mass of the virtual 
photon 7*, E its energy and p = |p| its momentum. Be- 
cause of M 2 = E 2 — p 2 > 4m* the rate is suppressed 
exponentially for temperatures below 2m Al . 

In order to obtain the spectrum from this formula one 
has to integrate over the space-time volume, taking into 
account the space-time evolution by using, for example, 
a hydrodynamical model. The total muon yield then fol- 
lows from integrating the spectrum over the energy and 
momentum of the virtual photon. 

At temperatures above 10 MeV also hadron production 
becomes important, in particular pion production |28j . 
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